High resolution SIMS depth profiling of nanolayers 
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We report results of high-resolution TOF SIMS (time of flight secondary ion mass spectrometry) 
depth profiling experiments on a nanolayered structure, a stack of 16 alternating MgO and ZnO ~5.5 
nm layers grown on a Si substrate by atomic layer deposition. The measurements were performed 
using a newly developed approach implementing a low energy direct current normally incident Ar + 
ion beam for sample material removal by sputtering (250 eV and 500 eV energy), in combination 
with a pulsed 5 keV Ar + ion beam at 60° incidence for TOF SIMS analysis. By this optimized 
arrangement, a noticeably improved version of known dual-beam (DB) approach to TOF SIMS 
depth profiling is introduced, which can be called gentleDH. 

We apply the mixing-roughness-information model to detailed analysis of experimental results. 
It reveals that the gentleDB approach allows ultimate depth resolution by confining the ion beam 
mixing length to about 2 monolayers. This corresponds to the escape depth of secondary ions, the 
fundamental depth resolution limitation in SIMS. Other parameters deduced from the measured 
depth profiles indicate that a single layer thickness equals to 6 nm so that "flat" layer thickness d 
is of 3 nm and interfacial roughness a is of 1.5 nm thus yielding d+2-a=6 nm. In essence, we have 
demonstrated that the gentleD~B TOF SIMS depth profiling with noble gas ion beams is capable of 
revealing structural features of a stack of nanolayers, resolving its original surface and estimating 
the roughness of interlayer interfaces, which is difficult to obtain by traditional approaches. 
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INTRODUCTION 



Secondary ion mass spectrometry, SIMS, is a well- 
established analytical method in surface, materials and 
chemical sciences. Besides conducting routine elemen- 
tal, molecular, and isotopic analyses of solid surfaces and 
thin films, creative applications of this method help in 
finding solutions to a variety of non-standard materials 
science problems [U-El]. Sputter depth profiling analysis 
provides information on the variation of sample compo- 
sition with depth below the initial surface. In such ap- 
plications, the SIMS method is in a class of its own due 
to the unique combination of high sensitivity and depth 
resolution. While for other surface sensitive analytical 
techniques, such as Auger Electron or X-ray Photoelec- 
tron Spectroscopies (AES or XPS), the inherent escape 
depth is >1 nm, for SIMS this characteristic parameter 
on the order of 1 to 2 atomic layers, which can be consid- 
ered as a theoretical limit of its depth resolution 0, 0] ■ 
In real life, however, the depth resolution of SIMS is lim- 
ited by an native surface roughness of the sample (if any) 
and by the fundamental effect called ion beam mixing [8| . 
Mixing is characterized by the ion beam mixing length w, 
which is a figure-of- merit of a SIMS instrument. For all 
depth profiling analyses, the former limiting factor is out 
of the analyst's control. Fortunately, one can minimize 
the latter effect by optimizing analytical approaches. Of- 
ten, new combinations of techniques for sample material 
erosion and composition probing are sought after. As 
an alternative to sputtering, sample material can be re- 
moved for compositional probing by laser enhanced field 



desorption 0, [l(| , as implemented in modern atom probe 
tomography instruments. In the case of SIMS depth pro- 
filing analysis, all such approaches, in general, aim at de- 
creasing the energy deposited into the sample by incident 
projectiles. This can be achieved by either directly de- 
creasing the ion beam energy, as usually done for conven- 
tional primary ions (Cs + , noble gas ions, and as well) 
111 Il2j | , or implicitly, by having more constituent atoms 
in the primary molecular/cluster ions (single charged SFg 
[3, G 60 molecules/Au and Bi [l|, or Ar llfj clus- 
ters) and dividing this energy between the constituents 
so that each atom has much lower impact energy (~l/60 
of the nominal beam energy per single constituent atom 
in the case of Ceo)- 

The dual-beam (DB) depth profiling [T^] was intro- 
duced in time of flight (TOF) SIMS instruments to de- 
couple the depth resolution of the direct current (dc) 
sputtering ion beam from the lateral resolution of the 
pulsed analysis ion beam. In the DB approach, the sput- 
tering and analysis cycles are independent, and the depth 
resolution is controlled only by the sputtering beam pa- 
rameters 18j. In order to fully realize the potential of 



the DB sputter depth profiling approach, we further im- 
proved it by implementing a normally incident sputter- 
ing ion beam, which is extracted from the ion source at 
nominal energy of a few keV and delivered with the same 
energy into the target region, where the target potential 
can decelerate it so that the impact energy does not ex- 
ceed a few hundred eV [l9j . This permits direct adjust- 
ment of the sputtering beam energy (down to the sput- 
tering threshold) without losses of the ion beam current 
and degradation of the beam spot quality. An important 
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benefit of the normal incidence sputtering is the signif- 
icant reduction in roughening of the surface induced by 
the ion beam [20j . Moreover, the lower the impact en- 
ergy, the lower the ion mixing effects are. Thus, the low 
energy normally incident sputtering ion beam maximizes 
the benefits of the DB sputter depth profiling and al- 
lows us to create a versatile instrument, suitable both 
for ultra-trace surface analysis 3, 21 1 and applications 
in the materials science applications, and thus address a 
current problem on successful mass application of nan- 
odevices discussed in Ref.[22J. To resolve this problem, 
development of advanced and powerful characterization 
tools, combining multiple points of view (i.e. probes) on 
a sample, is highly needed to achieve faster, facile and 
standardized characterization and so to reap the benefits 
of nanotechnology [22] ] . 

We apply here gentleUQ approach to better under- 
stand chemistry and structure of ultrathin films produced 
by atomic layer deposition (ALD). The ALD technique 
is widely used for a layer-by-layer materials synthesis 
and has great application potential in many areas due to 
its ability to coat high-aspect-ratio substrates by confor- 
mal layered structures with precisely predefined composi- 
tional profiles [23, 24 1 . One of such important application 
areas is functionalization of material surfaces for use in 
novel detectors and sensors. Our effort here is the con- 
tribution to the Large Area Picosecond Photo-Detector 
(http://psec.uchicago.edu/) collaboration, which is fo- 
cused on fast particle detectors with large areas function- 
alized by ALD. It is known that thin films of metal oxides 
can serve as very efficient emitters of secondary electrons 
and be used for improving detector performance. More- 
over, it has been recently shown that such ALD-grown 
films of several nm thickness enable independent tailor- 
ing of the electrical resistance and the secondary elec- 
tron emission yield 25] . A functionalization of a detec- 
tor based on a large area (tens square centimeters) mi- 
crochannel plate means that surfaces of the vast amount 
of microchannel plate pores must be coated by thin films 
with predefined characteristics. This task is the perfect 
match to the unique capabilities of ALD. At the same 
time, understanding and controlling the coating unifor- 
mity (thickness fluctuation etc.) and the condition of in- 
terfaces between resistive and emitting layers as well as 
between those and the substrate are extremely important 
for improving the materials' performance and function- 
ality. 

With this goal in mind, we characterize in this work a 
stack of 16 alternating ~5.5 nm MgO and ZnO layers (8 
of each) using gentleDB SIMS depth profiling approach 
with normally incident 250 eV and 500 eV Ar + ion beams 
for sputtering and a pulsed 5 keV Ar + beam (pointed at 
60° from the target normal) for analysis. The chosen 
stack of 8 pairs MgO/ZnO, being one of the proposed 
functionalization systems, serves two purposes. The first 
one is to create a planar (unfolded) model of an actual 



microchannel plate pore surface coating with complex ge- 
ometry coating in order to reveal the layer-to-layer repro- 
ducibility of the ALD process starting from the substrate, 
along with the layers uniformity and the conditions of in- 
terlayer interfaces. The other purpose is to test the depth 
profiling capabilities of the gentleDB approach. Besides, 
multilayer nanostructures and devices are at the forefront 
of materials science. Giant magnetoresistance multilayer 
films [2^] or superlattices of transition metal oxides [27j 
are striking examples of this. Desired performance of 
such structures strongly depends on the interfacial rough- 
ness, interdiffusion between layers, layer-to-layer consis- 
tency, and layer conformality as well. Thus, we aim to 
demonstrate in this work the unique synergy of the com- 
bination of the low energy sputtering with inert gas ions 
and the normal incidence angle, which makes the gen- 
tleUQ depth profiling approach superior to many others, 
especially when applied to characterization of nanolay- 
ers. In this case, the capability of sputter depth profiling 
without alteration of the sample chemical composition 
and modification of its surface and interface morphol- 
ogy does minimize instrumental artifacts and experimen- 
tal/procedural assumptions and thus helps to elucidate 
the relationship between these parameters and advance 
materials synthesis approaches. We will support these 
statements by demonstrating accurate depth profiling on 
the nanometer scale and determining key characteriza- 
tion parameters such as the ion beam mixing length, and 
the decay Xd and leading X up lengths, which effectively 
characterize the resolution of an instrument and can be 
compared with those of other approaches/setups in the 
field. The proposed characterization approach becomes 
even more important if one imagines that less-than-nm 
depth resolution is achieved laterally over the area of 
~1 mm 2 (in general, this area can be varied between 
^100x100 /^m 2 to several mm 2 ). 



SAMPLES 

A layered structure |5.5 nm MgO/5.5 nm ZnO|x8 was 
grown by the ALD technique on a Si substrate and char- 
acterized by x-ray diffraction and ellipsometry. The stan- 
dard calibration and chara cteriza tion procedures of ALD 
can be found in Refs.124 
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The structural data 
obtained from these measurements indicated that MgO 
layers are amorphous and ZnO layers are polycrystalline 
in the wurtzite phase. Although the surface roughness 
of the MgO/ZnO sample was not measured, Al 2 03/ZnO 
multilayers prepared under similar conditions had a sur- 
face roughness of 0.9-1.5 nm 24| 3(J. The roughness of 
the Si substrate was ~0.3 nm. The initial layer-to-layer 
mixing due to thermal diffusion during growth at T—A73 
K is expected to be extremely low. 
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EXPERIMENTAL 

The depth profiling of the layered structure described 
above was performed using a custom made SARISA 
(Surface Analysis by Resonance Ionization of Sputtered 
Atoms) instrument operated in SIMS mode (instead of 
its primary mode of secondary neutral mass spectrometry 
with laser post- ionization) []J, |3l|. The SARISA TOF 
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mass spectrometer was designed to operate in multiple 
modes using the same set of ion optics With post- 
ionization lasers switched off, the instrument operation 
corresponds to that of TOF SIMS with long primary ion 
pulses (~0.2-l /is) and delayed extraction of secondary 
ions. For the experiments described below, SARISA's 
ion optics either delivered the primary sputtering Ar + 
ion beam to the sample assuring its orthogonal impact 
onto sample surface and controlling its impact energy, or 
it extracted secondary ions generated by a pulsed analyt- 
ical Ar + ion beam and performed their TOF MS analysis 
[l9j . Fast switching between these regimes was done elec- 
tronically by changing the potentials of the optics elec- 
trodes. 

The low energy ion beam is formed by injecting a keV 
ion beam produced by a VG EX05 gun into the TOF 
spectrometer ion optics, by delivering and focusing the 
beam at normal incidence to the sample and by deceler- 
ating the ions to the desired low impact energy by biasing 
the target with an appropriate voltage (see Fig.l). The 
beam defocusing, which accompanies this deceleration, is 
compensated by changing the voltage of the outlet elec- 
trostatic lens called Lensl [19|. The impact energy of 
this milling beam was set in this study to either 250 or 
500 eV, and the ion current was ~1 fiA. This beam was 
raster scanned over the sample surface by engaging an 
octupole deflector present in the TOF system ion path 
(the Shaping Octupole in Fig. 1), so that a square crater 
of ~ 1.5x1.5 mm 2 (shown in Fig.l as a blue square) was 
eroded. The exact size of the crater depends on the decel- 
eration potential of the target, which additionally deflects 
the beam when it is raster scanned. For example, the ra- 
tio of the crater sides between £(,=250 eV and £&=3 keV 
equals 1.3; this value was proved both by crater imaging 
and by SIMION 3D© 32] simulation. 

The analytical ion beam, used to probe the sample 
at various depths, comes from an Atomika WF421 gun 
pointed at 60° incidence with respect to the target nor- 
mal (Fig.l). This ion beam can be independently tuned, 
raster scanned, and pulsed. The primary Ar + ion energy 
was always 5 keV, and the pulse duration used in the 
analysis in all the experiments was 200 ns, since mass 
resolution under these conditions was sufficient. The 
raster scan size of this ion beam in all the experiments de- 
scribed here was set to a 500x500 /im 2 (the green square 
in Fig.l). During the TOF MS analysis cycle, the de- 
celeration potential on the target is switched off. In the 




FIG. 1: (a) Diagram of the dual-beam system. The primary 
Ar + beam comes from the VG low energy column and is de- 
flected into the front TOF column by the Bending System 
optics 19]. The ion beam is tightly focused with the front 
end Lensl, impinging at normal incidence with respect to the 
sample surface. The beam is raster scanned by the Shaping 
Octupole. The separate Atomika ion gun is pointed at the 
target at 60° from its normal and equipped with a pulsing 
and raster scanning plates. The electron gun images surface 
topography, (b) White light interferometry image of a high- 
energy Atomika Raster crater positioned in the center of a 
low-energy VG Raster crater. The Atomika Raster crater was 
dc sputtered just after the completion of the gentleDB SIMS 
measurements in order to reveal the place where the analysis 
was done, since the pulsed analytical beam does not produce 
any visible bowl under conditions satisfying a <C 1. 



gentleDH depth profiling regime, the Atomika ion beam 
was only used for analysis. By choosing an appropriate 
aperture, the primary ion current was set to 30 nA, in 
order to attenuate the secondary ion signal so that the 
TOF MS detector is not overwhelmed. 

For the two ion beams used in these experiments, pa- 
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In order to perform the gentleDH depth profiling, we 
fit the analysis beam raster into the center of the low 
energy sputtering beam raster. To do this pre-alignment, 
we used a Faraday cup for ion beam collection and took 
advantage of available crater imaging capabilities. Crater 
imaging was done in situ by an optical Schwarzschild- 
type microscope [l9j and the SE M fToL [33| or by ex situ 
optical white light interferometry [34( (Fig. lb). 



RESULTS AND DISCUSSION 

Figure 2 demonstrates normalized depth profiles ob- 
tained by the gentleDB approach at 250 and 500 eV sput- 
tering beam energies. Both the 24 Mg and 64 Zn depth pro- 
files are very similar and feature high contrast and high 
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genlleDB profiles at Ar sputtering beam 

500 eV ™Mg | <> 250 eV M Mg 

500eV S1 Zn|— 250eV B,1 Zn 

500eV 2S S« |— o— 250eV 28 5i 

model ALD structure Tor " Mg profile 



Table I. Leading and decay lengths for 24 Mg, 64 Zn, and 
28 Si extracted from gentleDB depth profiles. 




c 

3 
■n 

[3 

co 

6 



-10 10 20 30 40 50 60 70 80 90 100 110 

Depth (nm) 

FIG. 2: SIMS depth profiles of the MgO/ZnO layers on a 
Si substrate in the gentleDH approach obtained by a 250 eV 
and 500 eV sputtering Ar + ion beam combined with a 5 keV 
analysis Ar + ion beam. The profile of the first MgO layer is 
poorly resolved because of surface contamination with Mg. 



stability of the signal all the way through the multilayer 
structure. The model structure of the 24 Mg distribution 
based on ellipsometry measurements is shown in Fig. 2 as 
a black dotted line. The end of the last layer depth profile 
is marked by a 28 Si signal spike from the substrate. The 
variation in the 28 Si SIMS signal is due to variation in 
the Si secondary ion yield 7 from SiC>2 , and its following 
drop and stabilization are observed at the 100-fold lower 
level with respect to the spike magnitude, which coincides 
with the known ratio 7 f,'° 2 ^100. The first peak in the 
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24 Mg profile has an unsymmetrical shape (cotangent-like) 
starting at the surface from the highest SIMS intensi- 
ties transforming right into trailing edge, while the 64 Zn 
profile starts approximately at the zero level; both pro- 
files then have a steady sine-like shape. Table I shows 
the leading X up and decay A^ lengths extracted from the 
depth profiles. The definition for the leading or decay 
length is the inverse of the first order coefficient for the 
linear approximation of the leading or trailing edges of 
a peak represented in semi-log scale, £n(Intensity) vs. 
depth [35j. Here, both X up and Xd are ~5 nm for 24 Mg, 
~2 nm for 64 Zn, and A„ p =0.4 nm for 28 Si. The mea- 
sured leading and decay lengths were on the same order 
or better than published to-date corresponding data ob- 
tained with commercial state-of-the-art SIMS instrumen- 
tation using cluster/molecular or low energy O J /Cs + ion 
sources, or with laser-assisted atom probe tomography 

a Ear 
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As mentioned above, the first 24 Mg peak of the stack 
was not resolved from the surface constituent, which caps 
the structure. In order to resolve in the depth profile 
the 24 Mg peak corresponding to the topmost ALD layer 
from the surface contamination, we decreased the time 
of the individual sputtering cycles and proportionally in- 
creased their number. Figure 3 shows the depth profile 
of the two top layers performed with more data points 
per depth unit to separate the 24 mass from the topmost 
layer 24 Mg peak. The green line here corresponds to the 
depth profile with poorer resolved Mg profile, shown in 
Fig. 2. This is to emphasize the capacity in the resolving 
power tuning of the gentleDH approach. The distance 
between subsequent depth points in Fig. 3 is ~0.2 nm, 
which precisely correlates with half of the average lattice 
parameter of 3.9 A for wurtzite ZnO [HI and half of the 
cubical unit cell size of crystalline MgO with a lattice 
parameter of 4.2 A [39]. 

The ratios of 25 Mg/ 24 Mg and 26 Mg/ 24 Mg proved that 
the mass 24 peak appears to be Mg with correct terres- 
trial isotopic ratio values. Such stoichiometry enrichment 
may be due to either natural surface contamination oc- 
curring after the sample is taken out of the ALD reactor, 
or due to some ALD specific effect manifesting itself at 
finalizing growth steps. At the same time, the high depth 
resolution profiling allows us to detect a "dip" between 
the surface 24 Mg and the topmost magnesia layer so that 
there is enough prior information to precisely reconstruct 
the contaminant and MgO layer profiles, assuming the 
contamination is natural. These two functions are dis- 
played in Fig. 3: cyan dotted line corresponds to MgO, 
while black dotted one is the surface Mg cap. 

These experimental results clearly demonstrate re- 
markable resolving power of the gentleDH approach to 
depth profiling. In this regard, keeping in mind that 
normal incidence should not affect pronouncedly native 
sample roughness [20j ] . let us qualify what the definition 
"gentle" means. For this, we apply the MRI (mixing- 
roughness-information) model [40j, |4l| to simulate the 
64 Zn peak in the 500 eV high resolution gentleDB depth 
profiles using the model parameters w (cascade mixing 
length), er (r.m.s., root mean square, roughness), and A 
(information depth, or escape depth, ~l-2 monolayer in 
the case of SIMS). 

A summary of the MRI model for the SIMS experiment 
is shown below. Note that the preferential sputtering 
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gentle~D~B profiles of high resolution 
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FIG. 3: High-resolution depth profiles of MgO/ZnO top layers 
at 500 eV sputtering Ar + beam and 5 keV analytical Ar + 
beam showing the surface Mg contamination being resolved 
from the top MgO layer peak. The green solid line is the 
starting part of the profiles copied from Fig.2. Cyan and black 
dotted lines are the Mg signals reconstructed from 24 Mg high 
resolution depth profile (open squares + dotted line in cyan) 
curve which correspond to 24 Mg in ALD magnesia and some 
sort of contaminant, respectively. 



effect is not taken into account here, since there is no 
such an effect for MgO and ZnO. As a starting point, the 
ion mixing is modeled first by the differential equation (I) 
for an ideal, flat layer of thickness d, with the boundary 
conditions (2), (4), and (5) 

cLCm{z) 1 „ , N 1 „ , . 
, + - • C M (z) = - ■ C act (z + to), (f) 
dz w w 



CZ(z)=A-(l-exp(-±+a) 
Cm(z)\ Zi -w = 



(2) 



Hence, 



A = 1 by default . (3) 

C^(z) = l-exp(-^±^) 

Cact(z + w)\ z > Z2 _ w =0. (4) 



Cm(z)\z 2 — w — CjJ[ (z)\ Z2 — w 



(5) 



Hence, 



B = f - exp(-&^f±) 

C d M {z) = (l-exp{-^))-exp{- 



(6) 



where Cm{z) and C ac t{z + w) are apparent (at any depth 
z) and actual (at w deeper relative to z) concentrations 
of an element, respectively; w is the ion mixing length; 
Z\ and z 2 are the flat boundaries of the arbitrary layer of 
interest, so that z 2 — Z\ = d equals the layer thickness, 
d. C^{z) and Cf { (z) are the leading and trailing edges 
of the depth profile peak, respectively. Thus, due to the 
cascade mixing, the solution yields 




z — zi +'UJ 



),z € [zi - w;z 2 



w 



expi-^ur^)) ■ exp(- 



Z — Z2+W 



),z> z 2 — w 



(7) 



If the surface roughness (either inherent, or induced, or 
both) is taken into account, then the profile additionally 
broadens as 



Cmr(z) 



1 



V2- 



C M {t) ■ exp(-%-^)dt, (8) 
2 • G 



where the mixing is convolved with the Gaussian func- 
tion that takes the r.m.s. roughness a as the standard 
deviation. 

Finally, when both the Cmr(z) parameter and the ef- 
fective depth contributing to the signal are considered, 
the depth profile shape can then be expressed as 



C M ri{z) 



Cmr(x) ■ exp(- 



X 



-)dx, (9) 



Z — Z2+W > 



where A is the information depth or ion escape depth, by 
conventional definition. Generally speaking, in the case 
of SIMS, it has negligible contribution in the elemental 
peak dependence on the depth Cmr(z), defined only by 
the cascade mixing and roughness. 

We found that the best matching between the mea- 
sured and modeled profiles occurs when the MM param- 
eters w and a were set as follows: w=0A nm, er=f.5 
nm, and the parameter A was fixed at value 0.2 nm cor- 
responding to 1 monolayer. The nominal thickness (flat 
thickness) of the layer, d, had to be of 3 nm to get the best 
peak fit in Fig. 4. The thickness d is not a free parameter 
of the MRI model. Nevertheless, being unknown a priori, 
it had to be varied in order to define the boundary con- 
ditions (2), (4), and (5). Then the effective thickness of a 
single ZnO layer, i.e., 2-a + d, equals to 6 nm. This value 
obtained by MRI does agree well with the ones estimated 
for a single ALD layer by ellipsometry and quartz crys- 
tal microbalance technique except the real material layer 
is not ideally flat, it can be represented by a structure 
made of flat inner layer squeezed between comparable in 
thickness additions on both sides representing wavy inter- 
faces. Thus, under gentleDB conditions, the value <r=1.5 
nm calculated for the first 64 Zn peak can be accepted as 
the native roughness of the sample inner interfaces (see 
Fig. 4). In agreement with the MRI model results, this 
roughness causes symmetrical peak broadening in sputter 
depth profiling experiments. 
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FIG. 4: MRI model fittings of the first 64 Zn depth profile 
peak: 500 eV, high resolution gentleDB profile - red circle 
data points vs. the MRI simulated curve in dark red. The 
inset is the same plot as the main figure, but in semi-log scale. 



It is interesting, that some preceding works, e.g. 



Ref . [42| , concluded that the attenuation of the signal in- 
tensity of each subsequent layer in laminate/multilayer 
structures, by approximately sin(ir) • exp(-x) law, oc- 
curs only due to the native interface roughness of layered 
structures. However, presented here results of the gen- 
tleDH depth profiling obtained on a sample of exceed- 
ingly disturbed grown interfaces suggest that the native 
roughness is not the dominant factor in the degradation 
and smearing of depth profiles in SIMS. To prove this, we 
applied the conventional single beam TOF SIMS depth 
profiling approach to the same sample, using 5 keV Ar + 
ions with 60° incidence angle (generated by the Atomika 
WF421 ion gun) - for both ion milling and TOF MS anal- 
ysis (results to be published). The comparison between 
single beam and gentleDH approaches revealed dramatic 
differences in results: the ion beam mixing length w in- 
creased from 0.4 nm to 3 to 4 nm (43|, which is on the 
order of the entire ALD layer thickness. Moreover, we 
observed with single beam approach the periodical sine 
signal with exponential attenuation mentioned above as 
well as the smearing of the peak shapes for deeper lay- 
ers. To summarize, the ion beam mixing appears to be 
the most important phenomenon to account for in depth 
profiling experiments on nanolayered structures. Reduc- 
ing ion mixing length w by a proper choice of the an- 
alytical procedure can produce depth profiles such as 
the ones shown in Fig. 2 that reveal much more infor- 
mation about the samples than about the measurement 
artifacts / distortions. 



CONCLUSION 

To conclude, we have outlined here our approach to 
solve the nanomaterial characterization bottleneck prob- 
lem [22[ in the materials science by application of TOF 



SIMS technique. Namely, we described our gentle vari- 
ant of the dual-beam sputter depth profiling (gentleDB) 
applied to characterization of nanolayered materials and 
demonstrated sub-nanomcter depth resolution of the 
method. This approach combines normally incident low 
energy (down to sputtering threshold) Ar + ion beam for 
sputtering with several keV pulsed Ar + ion beam (60° 
incidence angle) for analysis. The shallow penetration 
depth of these sputtering primary ions results in an ultra- 
short range of a crystal lattice disturbance of about 0.4 
nm thus enabling measurements of structural parameters 
(flat thickness and roughness) of a stack of 16 alternat- 
ing MgO and ZnO nanolayers (8 of each, with thickness 
of ^5.5 nm) grown by ALD on a Si substrate. The in- 
terfacial roughness as small as 1.5 nm estimated by the 
gentleDH sputter depth profiling does agree well with 
measurements obtained by other techniques: the surface 
roughness of individual ZnO layers (~0.5 nm) determined 
by atomic force microscopy [30( and the roughness of the 
multilayer structure (~0.9-1.5 nm) as determined by x- 
ray reflectivity and atomic force microsco py f or compa- 
rable ALD-grown Al 2 3 /ZnO multilayers [H ■ 
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